We used the fPrioritizrRun() function to streamline the runs. fPrioritizrRun() arranges all the data and makes sure that they’re in the right format, then runs problem() and solve(). This function also saves the summaries each plan (cost, % area protected and number of representative features that achieved 30% or more protection).
The relevant scripts for this step is found in ~/scripts/ with the prefix 10_Prioritizr.
# the function is in:
source("scrpits/10a_Prioritizr_fPrioritizr.R")
# to run the function with climate-smart data
source("scripts/10h_Prioritizr_fPrioritizrIUCNSmart_NoProv.R")
# to run the function with climate-uninformed data
source("scripts/10i_Prioritizr_fPrioritizrIUCNUninformed_NoProv.R")
# to create no-regret closures with the climate-smart plans
To create no-regret plans, we used thre function fCreateNoRegret() which intersected the selected PUs in all three climate scenarios. The relevant scripts for this step is found in ~/scripts with the prefix 11_NoRegret.
# the function is in:
source("scripts/11_NoRegret_fCreateNoRegret.R")
# the runs are found in:
source("scripts/11e_NoRegret_fCreateNoRegretIUCNRuns_NoProv.R")
As mentioned in the previous section, we will only be showing results of Runs 2 (Max = 0.9, Min = 0.1), 5 (Max = 0.6, Min = 0.1) and 8 (Max = 0.3, Min = 0.1).
SSP1-2.6 run
SSP126_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP126.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP2-4.5 run
SSP245_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP245.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP5-8.5 run
SSP585_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP585.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
noregret_target90 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target90.rds")
uninformed_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_uninformed.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 1 x 12
## target_scenarios sum_area percent_area num_pu total_cost median_velocity
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Target90_summaryuninf… 1.68e7 19.8 6308 120539. 16827221.
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## # median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## # scenario <chr>
SSP1-2.6 run
SSP126_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP126.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP2-4.5 run
SSP245_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP245.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP5-8.5 run
SSP585_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP585.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
noregret_target60 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target60.rds")
uninformed_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_uninformed.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 1 x 12
## target_scenarios sum_area percent_area num_pu total_cost median_velocity
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Target60_summaryuninf… 11203920 13.2 4200 73545. 11203920
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## # median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## # scenario <chr>
SSP1-2.6 run
SSP126_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP126.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP2-4.5 run
SSP245_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP245.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP5-8.5 run
SSP585_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP585.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
noregret_target60 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target60.rds")
uninformed_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_uninformed.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 1 x 12
## target_scenarios sum_area percent_area num_pu total_cost median_velocity
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Target30_summaryuninf… 5601960 6.58 2100 32055. 5601960
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## # median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## # scenario <chr>